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The aim of this paper is twofold. First, we introduce a new class of lin- 
earizations, based on the generalization of a construction used in polynomial 
algebra to find the zeros of a system of (scalar) polynomial equations. We 
show that one specific linearization in this class, which is constructed nat- 
urally from the QR factorization of the matrix obtained by stacking the 
coefficients of A{x), has good conditioning and stability properties. 

Moreover, while analyzing this class, we introduce a general technique to 
derive new linearizations from existing ones. This technique generalizes some 
ad- hoc arguments used in dealing with the existing linearization classes, and 
can hopefully be used to derive a simpler and more general theory of lineariza- 
tions. This technique relates linearizations to pencil arithmetic, a technique 
used in solving matrix equations that allows to extend some algebraic oper- 
ations from matrix to matrix pencils. 

1 Introduction 

The most used technique for finding the eigenvalues of a matrix polynomial is converting 
it to a linear problem. Given a matrix polynomial A(x) := X]f=o ^d^'^ £ C"'"[3;], a degree- 
1 matrix polynomial (also called matrix pencil) L{x) G C'^"''^"[x] is a linearization if there 
are E{x),F{x) G €'^"•'^"[2;] such that 

L{x) = E{x) diag(^(x), I(rf_i)„)F(x), (1) 

and det £'(x), det -F(x) are nonzero constants. Linearizations have the same eigenvalues 
and Jordan structure as the starting matrix polynomial [T^. Several different methods 
to construct such a pencil exist; we recall here the more common ones. 
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Companion forms [13] The well-known companion matrix of a scalar polynomial gener- 
alizes easily to matrix polynomials. An example of this construction is Co — Cix, 
where 



-Ad 



In 




/„ 




(2) 



In literature several variants of this form appear: the polynomial may be reversed, 
so that Ad (instead of ^0) appears alone on a coefficient, and the Ai can appear 
in the last block column instead of the first, or in the top or bottom block row. 

Companion forms in different polynomial bases [IJ The companion form can be adapted 
using in its construction different (orthogonal) polynomial families instead of the 
usual monomial basis, yielding for instance pencils in the form Co — Cix, where 



Ci := 



-Ai + (3oAo 
-A2 + ao^o 
-A3 



-Ad- 



71/ 

ail I32I -fsl 
02! PsI 



Id-il 
ad-2l Pd-il 



has a tridiagonal pattern related to the associated three-term recurrence. 

Vector spaces of linearizations [TH] Another generalization of ^ is considering the vec- 
tor space L2 of pencils Lq — Lix that satisfy 



Idn 
Odxdn 



Ln 



Idn 



Li 



^0 

Ai 
Ad 



(3) 



for some v G C^. The operation on the left-hand side is called row-shifted sum in 
|18j . A second vector space Li is constructed with a similar relation, generalizing 
another companion form in which the coefficients of the polynomial appear on a 
row. Pencils in the intersection DL := Li n L2 have many interesting properties; 
for any regular matrix polynomial, most of them are linearizations, and there is a 
simple test to tell when and to characterize their condition number. 



Fiedler pencils [2j Define the matrices Fi := diag(I(j_i)„, Cj, In[d-i-i)] 
where 

'a T 
— T n 



for i = 1, 2, . . . , d. 
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For each permutation ct of (1, 2, . . . , d), the pencil Cq— nf=i is a hnearization; 

in particular, the two permutations a{i) = i and a{d+l — i) yield companion forms, 
so this construction can again be regarded a generalization, though in a different 
sense, of the companion form. All the pencils obtained in this way can be laid out 
explicitly as a simple combination of blocks Ai, I and zeros. Several additional 
generalizations of Fiedler pencils exist I21j. 

Not all linearizations preserve the Jordan structure of the eigenvalue oo [1'6\ I16j : to 
ensure that {strong linearization property), we have to require additionally that revL(rE) 
is a linearization for rev^(a;), where the operator rev is defined by 

d d 
i=0 1=0 

Similarly, not all linearizations extend seamlessly to singular polynomials. Usually, the 
minimal indices of the singular part of the matrix polynomials are modified; therefore, in 
order to reconstruct the singular structure of A{x), we need to analyze how it is related 
to that of L(x) [TO]. 

Other additional properties are often desirable in a linearization. For instance, when 
the starting polynomial is structured (i.e., symplectic, or even, or Hermitian), we would 
like our linearization to have the same property |17J. Moreover, in addition to the 
eigenvalues, it is useful to be able to recover the eigenvectors of A{x) with little additional 
work (preferably none at all) from those of L{x) [Ml El- Last but not least, the condition 
number of an eigenvalue in L{x) can be higher than the corresponding one in A{x) |15j . 
leading to instability in the computation; it would be good to be able to control this 
condition number growth. 

The aim of this paper is twofold. First, we introduce a new class of linearizations, based 
on the generalization of a construction used in polynomial algebra to find the zeros of a 
system of (scalar) polynomial equations. We show that one specific linearization in this 
class, which is constructed naturally from the QR factorization of the matrix obtained 
by stacking the coefficients of A{x), has good conditioning and stability properties. 

Moreover, while analyzing this class, we introduce a general technique to derive new 
linearizations from existing ones. This technique generalizes some ad-hoc arguments used 
in dealing with the existing linearization classes, and can hopefully be used to derive a 
simpler and more general theory of linearizations. This technique relates linearizations 
to pencil arithmetic |5], a technique used in solving matrix equations that allows to 
extend some algebraic operations from matrix to matrix pencils. 

2 Basics 

For each C- vector space M (and in particular we are interested in M = C" or M = C"'") 
and for each degree d G N we define 

Md[A, ^] := {adX'^ + a^^iA^" V + • • • + ao/ : Oi e M}, 
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the space of degree-d homogeneous matrix polynomials in two unknowns with coefficients 
in M. We identify a single-variable matrix polynomials) A{x) of degree d with its 
projectivization A{X/ p)^'^. This makes it easier to deal with the eigenvalue x = oo (which 
becomes A = = 0); on the other hand, a purely projective theory of linearizations 
seems to be elusive, as a projective analogue of ([T]) and of the strong linearization 
property does not exist, up to our knowledge. 

Two matrix polynomials P{x),Q{x) G C"'"[2;], are called equivalent if there exist 
E{x),F{x) G C"'"[2;] such that E{x)P{x)F{x) = Q{x) and deiE{x), detP(x) are 
nonzero constant. They are called strongly equivalent if one can find constant matri- 
ces E{x) = E G C"'", F{x) = F £ C"." realizing the equivalence. Notice that both 
are equivalence relations. Any pencil equivalent to a linearization is a linearization; any 
pencil strongly equivalent to a strong linearization is a strong linearization. 

A matrix polynomial is called regular, or nonsingular, if det^(x) is not identically 
zero. Nonsingularity of a matrix polynomial A{x) = X]x=o ^i^^ implies that 



have respectively full column and row rank; the viceversa does not hold, though |14j . 

Theorem 1. Let a matrix polynomial A{X,fi) := X]f=o^j'^V^ * ^ C)^'"!-^)^] given, 
and let Co,Ci be as in 1^. Then, 

1. the pencil C{X,fi) := ACi — fJ-Co, known as second companion form, is a strong 
linearization ofA{X,fi). 

2. Let X and y* be right and left eigenvectors of A{X, fi) with eigenvalue (A, /i), i.e., 
A{X, fi)x = and y*A{X,fi) = 0. Then, the vectors 



Ai 



and A 



Ao A 



1 



Ad 



Ad 



X 



X = 



(Ai + l^Ao)x 
{A2 + + l^Ao)x 



{Ad^i + l^Ad-2 + ■■■ + ^Ao)x 



X 



-{^A2 + ^A, + ... + ^Ad)x 



(- 



Ad-1 + ^Ad)x 




y* = /^ly* fi'^'^Xy* 



X' 
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are right and left eigenvectors of ACi — I-iCq, respectively. Notice that the first 
expression for x is defined for A / 0, the second for /x 7^ 0, and they coincide when 
they are both defined. 

Lemma 2 (Bases completion lemma). Let W G ^■rn,m+n yjj^^j^ jy^n j,^,^^ rank and A £ 
C™+"'" with full row rank such that WA = 0. Then, one can find V e C"'"+™, B G 

(^n+m,m ^^^^ ^^^^ 



w 



A B 



Proof. Let us complete the columns of A and W* to a basis of C"*"*"", i.e., choose 13 £ 



and V e 



71, m+n 



such that 



A B 



and 



V* W* 



are nonsingular. Define 



Rii R12 




'v 


R22 




w 



A B 



Rii, R22 are invertible, since the right-hand side is the product of two invertible matrices. 
Finally, define B and V from the relations 



A B 



A B 



I —R12R22 

i?22^ 



V 
W 



R 



-1 
11 





V 
W 



□ 



3 Polynomial rings and Moller-Stetter algorithm 

In commutative algebra, a polynomial ideal / = (/i, ■■■,//;) Q C[xi, . . . ,Xd\ is called 
zero- dimensional if the system 

...,Xd) = 
fk{xi, ...,Xd) = 

has only a finite number of solutions, or, equivalently [S| Finiteness Theorem, Section 
2.5], if the space C[xi, . . . ,Xd]/I is finite-dimensional. When this holds, we have the 
following result. 

Theorem 3 (Moller-Stetter). Let I be a zero- dimensional ideal, and f G C[xi, . . . , Xd]/L. 
Consider the multiplication map Mg : C[xi, . . . , Xci]/L — )■ C[xi, . . . , Xd]/L defined as f ^ 
fg. When a basis o/C[xi, . . . ,Xd\/I is chosen, this map is represented by a matrix. Rs 
eigenvalues (counted with multiplicity) are the values of g{xi, . . . ,X(i) at every solution 
(counted with multiplicity) xi,...,Xd of the polynomial equations. 

Solutions of the polynomial system can then be found by solving the associated eigen- 
problem [9j. Even though this stronger condition is usually not of interest in the com- 
mutative algebra applications, in the 1-variable case it is possible to prove that Mx is 
a linearization of the single polynomial equation /i = that generates the ideal (there 
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is a single one, as C[x] is a PID). In fact, when the monomial basis is chosen, the 
multiplication operator is the companion matrix of fi. 

We aim to generalize this result to polynomial matrices and to the projective setting 
in order to produce linearizations. Given a projective matrix polynomial A{X, fi) := 
X]f=o^i^V^~* £ ^d'^^l^jfAt ^6 define the analogous of the "generated ideal" (A) as 

{A) := {v £ Q[A,)u] : ^(A,^) = A{X, fj.)u for some u G C"}. 

Note that (A) is a vector space. When the polynomial is regular, it has dimension exactly 
n (proof: the map u i— )■ A{X, fi)u must have trivial kernel); in fact, to prove this we only 
need the weaker condition that the Ai have trivial common (right) kernel. 

A natural generalization of the multiplication map used above is considering the 
maps 

Ma :/(A,^)^A/(A,^), 
M^:/(A,^)^^/(A,^). 

Notice, however, that these operations increase the degree of the homogeneous polyno- 
mials they operate on, thus it is not immediate to see what we should choose as their 
domain and codomain to replace C[xi, . . . , x^]//. After some attempts, one ends up 
with the idea of considering the quotiented operators 

Ma,Mm:Q-i[A,H^Q[A,/.]/(^), 

Mx:v{X,fi)^ Xv{X,is) + iA), (4) 
: v{X,fi) ^ f^v{X,n) + (A), 

and the hope that fJ,Mx — XM^ is a linearization. Such a hope is well founded, as we see 
in the following theorem. 

Theorem 4. Let A{X, fi) = A^X'^ + Ad-iX'^~^ + . . . + Aq^i'^ be any matrix polynomial 
such that (A) has dimension n. Let W £ C'^"''-'^"^^)" be any matrix with full row rank 
such that WA""^ = 0. 



L The matrices 



Wo := W 



^•n,dn 
Idn 



Wi := W 



Idn 



(5) 



represent, using suitable bases, the quotiented multiplication operators M\ and 
respectively. 

2. if A{X, fi) is nonsingular, then W{X, fi) := /iVFo ~ XWi, or more in general i^Mx — 
XM^ (for any choice of the bases) is a strong linearization. 

We give here only a proof of Part [TJ and defer the one of Part [2] to the next section, 
as it is immediate once we set up the right machinery. 



6 



Proof (of Part[J\ only). For each degree A; G N, we choose 



k-l • 



,k-2\2. 



as a basis of C^[A, /u]. Then, it is easy to represent (A) and the non-quotiented multiph- 
cation operators 

Ma : f(A,/i) At;(A,^), 
: f(A,/i) 1-^ /if(A,/x); 



namely, 



(yl) =Imy4' 



col 



We now need a basis of the quotient space C^fA, /u]/(^). We know from elementary 
linear algebra that if we complete a basis of (A) to a basis of C2[A,/x] in any way with 
fn+i, fn+2, ■ ■ ■ , /((i+i)n) then the vectors + (A) form a basis of the quotient space. 
Lemma El gives us V e C"''^" and B G C''"'" such that 



V 
W 



in particular, the columns of B complete those of A'^°^ to a basis and therefore induce a 
basis of the quotient space. In particular, since 

v + {A) = + BW^ v + {A) = BWv + {A), 

then Wv gives the coordinates of v{X, fi) + (A) in the chosen quotient space basis. In 
particular. 



AM. - /xMa =1^ A 



Idn 



(6) 
□ 



4 Dual linearizations 



We introduce in this section a general technique to derive new linearizations from existing 
ones. 

Theorem 5 (Left duality). Let L{\,fi) := fiLo — XLi G Ci''^ [X, fi] be a nonsingular 
matrix pencil, and let M{X,fi) := /uMq — AMi G Cf^'^[A,^] be another matrix pencil of 
the same size such that 

1. MiLo = MoLi; 

2. M™'" has full row rank N. 
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Then, M(X,iJ,) is strongly equivalent to L{X,fi). In particular, if the latter is a (strong) 
linearization of a matrix polynomial, then the former is a (strong) linearization as well. 

Notice that the first condition can be rewritten in the more appeahng form 



^JL^^^ =0, 



J:-- 





-In 



In 


' col 



Proof. We first prove that We apply Lemma [2 to M™™ and JL , obtaining 

l2n- 











Bi 


Mo 


-Ml 




Lo 


Bo 



Let us choose (a, (3) £ C^, normahzed with \af + 1/3|^ = 1, such that L{a, /3) is nonsingu- 
lar; then, we have the following identity, where * denotes blocks that are not of interest 
here: 



Li Bi 
Lq Bq 



L{a,(3) * 




aiN 




* * 




^In 


CtlN 






( 




Vi ' 




= 




Mo 


-Ml 





* * 




* M(a,/3) 



al]^ PIn 
—PIn olIn 



Standard Schur complement theory shows that M(a, /3) is the inverse of the Schur com- 
plement of —L{a,P), and thus the former is nonsingular whenever the latter is so. 

We could probably extend this argument to prove that the Jordan structures of the 
two pencils coincide; however, we switch to a different one that reveals more clearly the 
form of the strong equivalence relation between the two pencils. Let (a, /?) be so that 
both M(a, P) and L(a, /3) are nonsingular; thanks to the first condition, the identity 



holds, and thus 



M(a, P)L{X, n) = M(A, /u)L(a, /?) 
M{a,P)'^M{\,n)L{a,P) = L{X,n). 



(7) 
□ 



If M(A, n) satisfies the hypotheses of Theorem^ we say that it is a left dual of L{X, fi). 
Notice that a left dual is uniquely determined up to left multiplication by a nonsingular 
scalar matrix; in particular, its right eigenvectors are well-defined. 

Corollary 6 (Right duality). The same results hold if we replace the two conditions in 
Theorem [3| with 

1. LoMi = LiMo; 

2. M'^"' has full column rank N. 
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In this case, we say that M(A, uniquely determined up to right multiphcation by 
a nonsingular scalar matrix, is the right dual of L(A,/u). 

With Theorem O in place, the proof of Part [2] of Theorem H] is a one-liner. 

Proof of Theorem^ Part\^ WiCq = WqCi, thus W{\,^) is the left dual of a compan- 
ion form! □ 



5 Constructing duals 



The relation MiLq = MqLi has been studied estensively in the context of pencil arith- 
metic and inverse- free matrix iterative algorithms |19t El [TTl |S] • Two main techniques 
exist for constructing Mq.Mi G C™'*" starting from Lq,Li G C"*'*^. 

QR factorization OH] Construct the QR factorization 







'Qu 


Q12 




'r 


Li 




Q21 


Q22 








and take Mi = (5*2, Mq = —Q\2- practice, a QRP factorization should be used, 
since A*^"^ being close-to-rank-deficient is a concern here. 

Example 7. Though formally it is a different problem, the same strategy works 
for selecting a good W to use in ([S]). Construct a QR factorization 





<3oo 


Qoi • 


■ Qod 




'R 




QlO 


Qii . 


■ Qld 









Qda 


Qd\ ■ 


■ Qdd_ 




_0_ 



(all the blocks are n x n) and take 



W 



Qoi 
Qu 



Q02 
Q12 



Qod 
Qld 



Qdl Qd2 



Q 



dd 



Enforcing an identity block \E, "2(J] Suppose that the identity matrix is a submatrix of 



that 



. Then, we can select a permutation matrix 11 G 



-<2m,2m 



and X G 



such 



n 



Then, the identity 



-X I 
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holds, and thus we can choose 

Ml -Mo 



-X I 



GeneraUzing sUghtly, if a m x m submatrix Y of 
we have 



n 

Li 



-1 



n 



n 



and thus 



Ml -M, 



-ZY- 



I 

ZY-^ 
I 



is known to be nonsingular, 



Y 



n- 



Example 8. The simplest nontrivial pencil in the Fiedler family [21 Example 2.5] 
is 



(8) 



We choose H as the block permutation (1, 2, 3, 4, 5, 6) i-^ (1, 3, 6, 4, 5, 2), so that the 
leading block equals the identity matrix. Then, the method yields 








^0 







/ 








A 


I 












-A2 












/ 







I 






Ml -Mo 



-Ao 10 
-I -Ai 10 
I A2 A3 



which is, up to sign changes, a companion form. Conversely, if one did not know 
about Fiedler pencils at all, they could derive ([HD with the same technique as the 
right dual of this companion form. 

6 Relationship with other known linearizations 

6.1 Vector spaces of linearizations 

A linearization L{X,fi) := ^uLq — XLi is a right dual of W{X,fi) if and only if 



= WoLi - WiLo = W\^ 
Since ker W = A^°^, this holds if and only if 



Idn 



Idn 



Idn 



Idn 



Ai 
Ad 



Z 



(9) 
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for some Z G C"'°'". This generalizes ([3]), and thus all pencils in L2 are contained in the 
right dual of M^(A, [i). Analogously, Li can be obtained as left dual of another companion 
form. 

In [18j, it is proved that any nonsingular L(A, /u) satisfying ([3]) is a strong linearization. 
Using Theorem [51 we can replace nonsingularity with a weaker condition. 

Corollary 9. Let A(A, /u) he a nonsingular matrix polynomial, and L(\, fi) be any matrix 
pencil satisfying Q. Then, L is a strong linearization of A{X, fx) if and only if L'^°'- has 
full column rank. 

6.2 Fiedler pencils 

We have already shown in Example [5] that pencil duality can be used to obtain naturally 
at least one Fiedler pencil starting from the companion form. This is no accident: 
the approach used in [2j to prove that the Fiedler pencils are linearizations can be 
reinterpreted in our setting as a sequence of duality operations that transforms a generic 
Fiedler pencil F{X,fi) into C{X,fi). Namely, one can recognize that Lemma 2.2] 
there is essentially a special case of our Theorem O In particular, the result proved here 
provides a correct proof of its part (b) , which is obtained in ^ using the false implication 
M{X, n) singular M™™ is rank-deficient. 



7 Eigenvector recovery 

It follows from ([7|) that for each (a, (3) such that A{a, (3) is nonsingular 

r ■.=y*W{a,P)-\ 
X :=C{a, f3)x 



(10) 



are left and right eigenvectors of W{X,fi). Different choices of {a, 13) in (jlOp yield the 
same x and y, up to normalization; this is proved by the following computations, which 
highlight the role of the expression /3A — a/i. 



{/3X - aij)y* 



AV] B{^Wo - aWi 
AVI BW (/3 



a 



Idn 

^nxdn 



A-'y*] (/(d+i)n - AV) [^/3 



^dn 



a 



Idn 

^nxdn 



a 



Idn 



and thus 



y 



j3X — an 



A^y 



B. 
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For the right eigenvalues, one can verify that in the case A 7^ 0, /i 7^ 
X =C(A, fi)x = 





(/3A — afi) — 



1 



=(^A - a/i)- 




1 

X 




1 



1^ 



1 




A 
1 



1 
1 



^1 



1 













A2 




Ad. 



(11) 



and when A = or /i = the version without vanishing denominators is vahd as weU. 



8 Conditioning 



In this section, we wish to estimate the conditioning of the hnearization Ty(A, //) proposed 
here, using the formula from [15J 



KH^(A,/U) 



\wif + \f,n\wo\ 



1/2 



m\ \\y\\ 



We estimate separately the factors in the formula. Using an argument similar to \\.b\ 
page 1020], we rewrite the denominator as 



y 



(pWi + \Wq)\^^ X =y* [fiCi + ACo) x = y* (jiDxP - XD^P 



(a,m) 



Normalize now so that |Ap + |/i|^ = 1, and notice that we are always in one of the two 



cases 



< 1, 

< 1, 



|A| < l/^l ^ 
l/"l<|AN 
In both cases, we may estimate in (|11|) 

\\x\\ < \j3X- aii\V2T{d) 

where 



< V2, 

< V2- 



A 



col 



r(d) 



1 1 ••• 1 

1 ■•. : 

■•• 1 

1 
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An elementary linear algebra argument yields 



T{d) 



2 sin 



2d+2 



2d 

TT 



(12) 



so T{d) grows linearly with d. 
An easier estimate holds for y: 



\y\\ < 



1 



I/3A — afd] 



|A| 



A, /J 



with 



|A|' \xf-'\^^\ 



Putting everything together, we obtain 



\B\ 



A| 



A 



col 



\B\ 



Fll \\y\ 



(A 



The bound gets simpler if we choose W with orthonormal columns as suggested in 
Example [3 in this case ||Wo|| ^ 1; ||l^i|| ^ 1) ll^ll ^ Ij a-nd thus 



1 col 



Fll \\y\ 



r (jx DxP - XD^P 



Apart from the moderate factor ^/2T{d), the only difference with the condition number 



of the polynomial eigenvalue problem is that |A| 



A, /J 



1/2 



|A| 



^col 



is replaced by 

l^ol 
l^il 

\Ah\ 



there. This reveals that the only problematic case is the one in which the two vectors 
appearing in this expression are nearly orthogonal. In the quadratic case, if proper 
scaling is performed, this can happen only if ||^i|| ^» ll^oll + ||^2|| and the eigenvalue is 
close to either or +oo. Incidentally, [15\ Theorem 4.3] shows that this is the problematic 
case for the BL linearizations as well. 

In fact, for W{X,^i) the conditioning of the linearization is only part of the story, 
since getting from the original data (coefficients of A{X, fi)) to W{X,n) requires a QRP 
factorization. Multiplying the conditioning of the QRP factorization by this condition 
number is likely to overestimate the error. Therefore, in the next section we report 
numerical results to obtain an practical assessment of the accuracy of the eigenvalues. 
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9 Numerical experiments 



Our primary goal in this section is comparing the quaHty of the eigenvalues computed 
from VF(A, n) with that of two common choices for unstructured problems. One is 
the companion form ([2]), which is used in Matlab's polyeig; the other is the strategy 
suggested in [TJ], namely, 

• if for an eigenvalue |x| > 1 (i.e., |A| > then use the DL linearization with 

V = ei 

• otherwise (|x| < 1, i.e., |A| < |/i|), use the BL linearization with v = e^- 

Their theoretical result proving a bound on the eigenvalue condition numbers is valid 
only in the case in which there are no or oo eigenvalues, i.e., Aq and are nonsingular; 
we restrict the comparison to this case for now. 

As stated before, comparing only the condition numbers of the linearizations would 
ignore the errors deriving from the QRP factorization; therefore, we chose to compare 
directly the forward errors on the eigenvalues. The same comparison is used in |15j . 

A set of reference values is obtained by using Matlab's variable precision arithmetic. 
Unfortunately, as of version R2010b, none of the customary methods for computing 
generalized eigenvalues (polyeig, qz and the two-argument version of eig) can deal 
with variable precision inputs; therefore, we have to settle for an inferior approach in 
computing the reference values. Namely, we compute them using eig(GO/Gl), where 
Go — xGi is the v = ei linearization. Both inversion and eigenvalue computation are 
performed with the default (for VPA) value of 32 significant digits. Incidentally, this is 
another reason to restrict our comparison to problems without infinite eigenvalues. 

Then, using the two-argument eig, we compute the eigenvalues of the pencil obtained 
by linearization with all the four linearizations under examination, and compare them 
to the reference values using the angular distance {d{v,w) =subspace( [v; 1] , [w; 1] )). 

The benchmark problems used are chosen among the regular quadratic eigenvalue 
problems {d = 2) in the collection NLEVP [6J, preprocessed with the Fan-Lin- Van 
Dooren scaling. As argued before, we restrict to those problems for which the condition 
numbers of and A2 are smaller than lO^'^; moreover, the slowness of variable precision 
arithmetic makes it impractical to deal with large problems, therefore we only consider 
the problems with n < 150. That leaves us with 17 problems, of which we have to further 
exclude problem signl as the computation of the reference eigenvalues did not converge 
due to a VPA error. We report the angular errors for the remaining test problems in 
the figures. 

The results are generally good, on par or better than the other alternatives. In prob- 
lems bicycle and metal_strip, the results are slightly inferior. On the other hand, 
in problems power_plant and wing the accuracy of the new method is higher, and mi- 
nor benefits on some eigenvalues can be identified in other problems. Moreover, notice 
that the strategy in [HMT] requires solving two separate eigenproblems, thus producing 
similar results with half of the computational work is an interesting result in itself. 
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We tested the new linearization VF(A, fi) on problems with several zero and/or infinite 
eigenvalues as well, without however obtaining good computational results. The reasons 
for this failure and the performance of this method in presence of zero and infinite 
eigenvalues are still under investigation, as it is not easy for us to identify the impact of 
the eigenvalue solver and of the various parts of the algorithm. 
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